Stability and entanglement in optical-atomic amplification of trapped atoms: the role 

of atomic collisions 
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Atomic collisions are included in an interacting system of optical fields and trapped atoms allowing 
field amplification. We study the effects of collisions on the system stability. Also a study of the 
degree of entanglement between atomic and optical fields is made. We found that, for an atomic 
field initially in a vacuum state and optical field in a coherent state, the degree of entanglement 
does not depend on the optical field intensity or phase. We show that in conditions of exponential 
instability the system presents at long times two distinct stationary degree of entanglement with 
collisions affecting only one of them. 

PACS numbers: 03.75.Be, 42.50.Ct, 42.50.Dv 
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I. INTRODUCTION 

Bose- Einstein condensation of trapped atomic gases [l|, 
H, [|[ has produced a fantastic advance in atom optics. 
Particularly, the interaction of condensates with single 
mode quantized lig ht fields has been a fascinating topic 

H 1, i, 0, 1, II Efl EL di 0, El HE M E3 , allowing 

for instance light and matter wave amplification [l2|, 
14 1, optical control of atomic statistical properties [la . 
16j, and potential applications in quantum information 
technology [l7| . 

It is known that the trap environment can modify the 
properties of ultra-cold atoms, such as its critical temper- 
ature fl8l . Il9j . Then it is expected that trap environment 
also influences the interaction between ultra-cold atoms 
and optical fields. In Ref. [2(J, the trap environment ef- 
fects on the Condensate Collective Atomic-Recoil Laser 
(CARL) [I!], [22| was considered by expanding the matter- 
wave field in the trap matter-wave modes. Such a situ- 
ation was called Cavity Atom Optics (CAO), in analogy 
with the Cavity Quantum Electrodynamics (CQE) where 
the spontaneous emission is modified by the presence of 
the cavity. The model obtained presents interesting prop- 
erties such as two regimes of exponential instability [20| 
and statistical properties that depend on intensity and 
phase of the optical field [23[. However, neither atomic 
collisions was taken into account nor a detailed study 
of the degree of entanglement between atomic and op- 
tical fields was done. Entanglement, one of the most 
notable characteristics of quantum mechanics (24|, has 
been object of intense study, both in systems involving 
light and matter [IE and in solids [13, HU , due to its 
important role for quantum- information processing [29| . 
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Therefore, inclusion of collisions and characterization of 
entanglement are important in order to turn the model 
more realistic and potentially useful for applications in 
the context of quantum information theory. 

This paper deals with the extension of the model 
of atom-optical parametric amplifier considered in Refs. 
[2fl |23| by including collisions between atoms. A study 
of the changes due collisions on the thresholds of expo- 
nential instability as well as in growth rate of the fields 
amplitudes is presented. Also a characterization of the 
degree of entanglement between atomic and optical fields 
is done. We show that in conditions of exponential in- 
stability the system presents at long times two distinct 
stationary degree of entanglement with collisions affect- 
ing only one of them. 

The article is organized as follows: In Section [TT] we de- 
rive the effective hamiltonian describing the system stud- 
ied. In Section IIIII we present an analysis of the system 
stability. In Section llVl we consider the atom-photon de- 
gree of entanglement in the regime of field amplification. 
Finally, in Section [V] we present the conclusion. 



II. MODEL 

We consider a Schrodinger field of bosonic two-level 
atoms, with transition frequency v, interacting via two- 
body collisions and coupled by electric-dipole interaction 
to two single-mode running wave optical fields of frequen- 
cies u>i and L02, treated as a quantum and a classical 
field, respectively. Both optical fields are assumed be- 
ing far off-resonant from any electronic transition. Thus, 
although the atom's internal state remains unchanged, 
the center-of-mass motion may change due the atomic 
recoil induced by two-photon virtual transitions. In the 
far off-resonance regime the excited state population is 
small, and therefore spontaneous emission as well as col- 
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lisions between excited state atoms, and between ground and the ground state atomic field evolves coherently un- 
state and excited states atoms, may be neglected. In this der the effective Hamiltonian 
regime the excited state can be adiabatically eliminated 
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where 

H = -^V 2 + V(r), (2) 
2m 

in is the atomic mass, V(r) is the trap potential, A = 
LU2 — v is the detuning between atoms and the classical 
optical field, g\ and g 2 are the optical fields coupling co- 
efficients, K = ki — k2 is the difference between their 
wavevectors, and S = lj\ — lo-i is the detuning between 
them. The operator Si is the photon annihilation oper- 
ator for the quantized optical field, taken in the frame 
rotating at the classical field frequency loi- The optical 
field treated classically is assumed to remain undepleted, 
so that d2 is simply a constant related to its intensity. 
Terms corresponding to the spatially independent light 
shift potential were neglected, so that the index of re- 
fraction of the atomic sample was assumed the same of 
the vacuum. Two-body collisions were included in the s- 
wave scattering limit by the second term inside brackets 
in the Hamiltonian (fT]), where 

U-^, (3) 
m 

and a is the s-wave scattering length, which depending 
on the repulsive or attractive character of the interaction 
can assume positive or negative values, respectively. We 
consider only positive values of scattering length, which 
are suitable for the creation of large condensates, and 
correspond to a situation more consistent with the ap- 
proximations in this paper. 

We assume that the atomic field is initially a Bose- 
Einstein condensate with mean number of condensed 
atoms N, and that this condensate is well described by 
a number state so its initial state is described by 

m = 0)) = -L=(cl) N |0) (4) 



where |0) is the vacuum state, and 

4 = J d 3 r ^ (r)* f (r) (5) 

is the creation operator for atoms in the condensate state 
¥>o( r )- Due to the presence of collisions the condensate 
wave function (po(r) satisfies the Gross-Pitaevskii equa- 
tion dJ[3i| 

(Ho + NU \Mr)\ 2 ) M*) =wW, (6) 

where fi is the chemical potential. 

Now we expand the atomic field operator in terms of 
the trap eigenmodes {ip n (r)} according to 

%(r) = tp (r)ca+64>(r), (7) 

where Sip(r) — J2^of n ( r )^ n crea t es particles in the 
excited trap modes, c„ is the annihilation operator for 
atoms in mode n, and J d 3 r ^(r)^J n (r) = S mn . We are 
interested in the linear regime, valid for interaction times 

so that J d 3 r S^Sr/j^ -C ^ejeo^. Then we can maintain 

only quadratic terms in the fields operators and invoke 
the undepleted approximation, which permits to substi- 
tute the condensate mode by a c-number evolving as 

c 0) « VNe-^. (8) 

Therefore, in the linear regime, by inserting expansion 
([7]) and taking into account the Gross-Pitaevskii equation 
(Eq. (|6|)), the Hamiltonian (fT]) reduces to 



H = TiSalai + hjd 3 r 6ft (H + 2NU) Sip + ^Ne 2i ^ J d 3 r pfSip 2 + h.c^j 
'-VNa\J d 3 r (e v V5e~ lK ' r <¥ + e** &ft e**" <p^ + h.c. 
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which in terms of the trap excited modes expansion be- 
comes 



III. STABILITY 



The solution of the linear system (|15[) can be written 



as 



+h (a + a f ) Xn ( 



(10) 



where 



fVjJr. 



J d 3 r<p* n (n + 2NU Mr)| 2 ) <p n , (11) 



is the collision modified energy of the nth trap level, 
Xn = VNA 0n \g x \ \g 2 \ \a 2 \ / |A| is the effective coupling 
constant between the condensate and the quantum opti- 
cal field, a — (gig^a^A/ \gi\ \g 2 \ \a 2 \ |A|)fii is the optical 
field operator multiplied by a phase factor that is related 
to the classical optical field phase, 



and 
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are the collision parameter and the element of matrix for 
the optical transition, respectively, both assumed being 
real numbers. The phase factors e ±v * were included in 
the operators c„ and c£. 

For simplicity, we assume that the overlap-integral in 
the collision parameter given by Eq. (| 1 2[) has a significant 
contribution for a given n = I = m and the matrix for 
optical transition Aq 71 is sharply peaked for n — m [37j . 
In this case we can neglect all excited trap modes except 
the m-mode in the Hamiltonian (fl"0|) . and the following 
effective Hamiltonian is obtained 



+ TlXm (a^cln + fit Cm + clnO- + ) ( 14 ) 



Terms like a^cj^ in Hamiltonian (fT4|) correspond to the 
generation of correlated atom-photon pair. These terms 
are analogous to the Nondegenerated Optical Paramet- 
ric Amplifier [HJ . Inclusion of collisions introduces terms 
like cj, 2 , which are responsible by creation of an atomic 
pair, and are analogous to the Degenerated Optical Para- 
metric Amplifier [321 ] . The Heisenberg equations of mo- 
tion for the field operators obtained from Hamiltonian 
(fT4|) result the following 4x4 linear system of equations 
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where the index m was dropped in order to simplify no- 
tation, and we introduced the dimensionless parameters 
t = oj m t, S = S/u> m , k = K m /uj m , and x = Xm/u m - 



3=1 



(16) 



where we defined x\ = c, x 2 = c\ £3 = a and X4 — for 
convenience, G t] (t) = ELi I U k \U~ 1 ]kje i " kt > [U]« * 
the ith component of the kth eigenvector of the matrix at 
RHS of Eq. (fT5|) . and ojk are the system eigenfrequencies. 

Stability analysis shows two regimes of exponential 
instability: (i) For S > and x 2 > 8(1 + k)/4 there are 
two purely real and two purely imaginary eigenfrequen- 
cies of the form {o>i = Q,uj 2 = — SI, £^3 = iT^uj^ = — iT}, 
where f2 and T are both real quantities. There is only 
one exponentially growing solution at the imaginary 
frequency W4 and the system is unstable. (ii) For 
S < and X 2 > (1 - k 2 - S 2 ) 2 /16 \S\ (1 - k) the 
eigenfrequencies are complex numbers of the form 
{u>i = Q + iT, uj 2 = -O. + iT, w 3 = Q, - iT, uj4 = -Q.- iT}. 
This case presents two exponentially growing solutions, 
W3 and u>4, which grow at the same rate T, but rotate at 
equal and opposite frequencies ±f2, producing a beating 
in the exponential growth of the fields intensities. 
Otherwise the eigenfrequencies are real and the system 
is stable. 

In Fig. [T] we plot in the 6 — x 2 plane the values of pa- 
rameters defining the threshold between stable and un- 
stable solutions. Points inside region I correspond to one 
exponentially growing solution, whereas for points inside 
region II there are two counter-rotating exponentially 
growing solutions. The full line indicates the threshold in 
the absence of collisions whereas dashed and dotted lines 
shows the change due collisions, which have the effect 
of reducing the unstable regions. In addition, for small 
X 2 the region II is centered around 5 = —y/l — K 2 and 
at asymptotically large detunings its threshold grows as 
X 2 > ~(5 3 /16 (1 - k). 

In order to illustrate the effects of collisions on the 
fields intensities 7j(t) = (^x\(t)xi(t)^, with % = 1,3, we 

consider that the atomic mode begins in a vacuum state 
|0) whereas the light field is initially in a coherent state 
|a). Then, with the help of Eq. lfT6|) , we found 

h(t) = \G l2 (t)\ 2 + \G l4 (t)\ 2 + \G i3 (t)a + G l4 (t)a*\ 2 . 

(17) 

Figs. [H-(a) and [21- (b) shows the logarithmic plot of the 
optical field intensity Is(t) as a function of time for pa- 
rameters lying in the regions I and II of instability, re- 
spectively. The intensity of the atomic mode has a sim- 
ilar behavior. We see that collisions reduce the rate of 
exponential growing and change the beating frequency of 
oscillation of the field intensity. The change in the den- 
sity distribution of atoms in the condensate mode due 
to collisions reduces the scattering of photons into the 
optical field, as well as the scattering of atoms by the 
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FIG. 1: Instability domains. Points inside region I correspond 
to one exponentially growing solution, whereas for points in- 
side region II there are two counter-rotating exponentially 
growing solutions. The full line define the threshold in the 
absence of collisions (k = 0.0). Change in the threshold due 
collisions correspond to the dashed [n = 0.4) and dotted lines 
(ft = 0.8). The parameters are dimensionless. 



optical field into the atomic mode, in opposition to field 
amplification. 



IV. ENTANGLEMENT 

Now we turn to analyze the atom-photon degree of en- 
tanglement for the exponential growing regimes of fields 
amplitudes. It is known that for a two-component sys- 
tem in a pure state the degree of entanglement can be 
characterized by the entropy or purity of one of the sys- 
tem components [34], [35[. Such entanglement measures 
require the calculation of the time-dependent quantum 
state for the system, which may not be an easy task. 
Instead, we consider a recently proposed entanglement 
coefficient in terms of cross-covariances of the fields op- 
erators defined by [36] 
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tors. Since we know the time dependent solutions for the 
fields amplitudes the parameter Y is easily calculated for 
the system considered in this paper. 

The coefficient Y satisfies the inequality < Y < 1, 
where at the maximum value of Y the system is maxi- 
mally entangled and was introduced taking into account 
that for two systems with operators A and B, if any of 
the cross-covariances with these operators differs from 
zero, then the modes are entangled. For the excited trap 
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FIG. 2: Logarithmic plot of the light field intensity as func- 
tion of time, (a) for parameters lying over the region I of the 
instability domains, (b) for parameters lying over the region 
II of the instability domains. Full line corresponds to absence 
of collisions (ft = 0.0) while dashed and doted lines corre- 
sponds to inclusion of collisions with ft = 0.4 and ft = 0.8, 
respectively. We set \ = 1 an< i ct = 2. The parameters are 
dimensionless. 



with the notation XiXj — (xiXj) — (xi) (xj) for the un- 
symetrized centered second-order moments of the opera- 



mode starting in a vacuum state and the light field in a 
coherent state, we obtain, 
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FIG. 3: Plot of the entanglement coefficient as a function of 
time in absence of collisions for several values of detuning S. 
We set x = !• The parameters are dimensionless. 



Eq. (fT5)) shows that the degree of entanglement do not 
depend on the optical field intensity or phase. Further- 
more, in the regime of exponential instability the param- 
eter Y attains at long times a stationary value thats is 
dependent on the signal of the detuning S. To see that, 
in Fig. [3] the parameter Y is plotted as a function of time 
for several values of detunings in the absence of collisions. 
We observe that for 5 > (one exponentially growing so- 
lution) the system attains at long times the maximum 
degree of entanglement, while for S < (two exponen- 
tially growing counter-rotating solutions) the long time 
degree of entanglement attains an oscillating stationary 
value far below the maximum one. 

The effects of collisions are presented in Figures 0J- (a) 
and HJ-(b) that show the entanglement coefficient as a 
function of time for values of parameters lying in regions 
of one and two exponential growing solutions, respec- 
tively and considering different values for the collision 
parameter k. We see in Fig. [4]- (a) that collisions do not 
affect at long times the degree of entanglement for 6 > 0, 
which always tend to the maximum value. However, for 
5 < (Fig. Hl-(b)) collisions strongly affects the long 
time degree of entanglement by changing the amplitude 
of oscillations of its stationary value. 



collisions reduce the growth rate of the fields amplitudes. 
For an atomic field initially in a vacuum state and optical 
field in a coherent state, we have verified that the degree 
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FIG. 4: Plot of the entanglement coefficient as a function of 
time, (a) for parameters lying over the region I of the instabil- 
ity domains, (b) for parameters lying over the region II of the 
instability domains. Full line corresponds to absence of colli- 
sions (k = 0.0) while dashed and doted lines corresponds to 
inclusion of collisions with k = 0.4 and k = 0.8, respectively. 
We set x = !• The parameters are dimensionless. 



of entanglement between atomic and optical fields does 
not depend on the optical field intensity or phase. Fur- 
thermore, in conditions of field amplification the degree 
of entanglement attains at long time a stationary value 
dependent on the regime of exponential instability, being 
maximum only in the case of one exponentially growing 
solution. Finally, atomic collisions only affect the long 
time degree of entanglement in the case of two exponen- 
tially growing solutions. 



V. CONCLUSION 
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